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A microscopic calculation of the equation of state for asymmetric nuclear 
matter is presented. We employ realistic nucleon-nucleon forces and operate 
within the Dirac-Brueckner-Hartree-Fock approach to nuclear matter. The 
focal point of this paper is a (momentum-space) G-matrix which properly 
accounts for the asymmetry between protons and neutrons. This will merge 
naturally into the development of an effective interaction suitable for appli- 
cations to asymmetric nuclei, which will be the object of extensive study in 
the future. 
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I. INTRODUCTION 

Nuclear matter is an idealized uniform infinite system of protons and neutrons under their 
mutual strong forces and without electromagnetic interactions. Symmetric nuclear matter 
(that is, equal densities of protons and neutrons) has been studied extensively. The so-called 
conventional approach to nuclear matter goes back to earlier works by Brueckner and others 
[1-6] and is known as the BHF (Brueckner-Hartree-Fock) theory. During the 1980 's, the 
Dirac-Brueckner-Hartree-Fock (DBHF) approach was developed [7-9]. The break-through 
came with the observation that the DBHF theory, unlike the conventional one, could de- 
scribe successfully the saturation properties of nuclear matter, that is, saturation energy 
and density of the equation of state (EOS). The DBHF method adopts realistic nucleon- 
nucleon (NN) interactions and contains features of the relativistic theory. It characterizes 
the nuclear mean field by strong, competing scalar and vector fields that together account 
for the binding of nucleons as well as the large spin-orbit splitting seen in nuclear states. The 
DBHF framework is a reliable, as well as feasible, microscopic method to describe effective 
interactions in the nuclear medium. 

Concerning asymmetric nuclear matter, systematic empirical investigations to determine 
its saturation properties have so far not been done. From the theoretical side, some older 
studies can be found in Refs. [10,11]. Interactions adjusted to fit properties of finite nuclei, 
such as those based on the non-relativistic Skyrme Hartree-Fock theory [12] or the relativistic 
mean field theory [13], have been used to extract phenomenological EOS. Generally, consid- 
erable model dependence is observed among predictions based on different EOS [14]. Varia- 
tional calculations of asymmetric matter have also been reported [15]. In Ref. [16], a Lorentz 
invariant functional of the baryon field operators is defined to project Dirac-Brueckner nu- 
clear matter results onto the meson-nucleon vertices of an effective density-dependent field 
theory. This is then applied to asymmetric matter and finite nuclei in Hartree calculations 
[17]. 
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We use realistic NN forces and operate within the DBHF framework. Similar calculations 
have been done in Ref. [18]. The purpose of our paper, however, is to produce a self-contained 
report, including all relevant formulas and nuclear matter results, and, most important, to 
examine the effect of the asymmetry on the in-medium NN interaction. This work will 
serve as a baseline for the future applications we plan to pursue. These include reaction 
studies using effective interactions based on asymmetric matter suitable for scattering on 
finite nuclei, and calculations of neutron radii and neutron skins. Concerning properties 
of astrophysical implications, we calculate the pressure in symmetric matter and neutron 
matter at high densities (up to about five times saturation density), and compare with recent 
experimental constraints obtained from analyses of nuclear collisions [19], with the extracted 
pressures being the highest recorded under laboratory conditions. 

Recently, much interest has developed around the study of highly asymmetric nuclei (in 
particular, extremely neutron-rich nuclei, and halo nuclei). If approved for construction, the 
Rare Isotope Accelerator (RIA) will map the limits of nuclear existence and allow the study 
of the unique nuclear systems which populate those boundaries. Thus, it is important and 
timely to develop microscopic effective interactions which can account for the asymmetry 
between proton and neutron densities. 

In the self-consistent approach, a calculation of nuclear matter properties yields, at the 
same time, a convenient parametrization of the density dependence in the form of nucleon 
effective masses. This information can then facilitate the calculation of the scattering matrix 
at positive energies, and finally of an effective interaction suitable, for example, for proton 
scattering on asymmetric nuclei. This interaction will be "isospin dependent" , in the sense 
of being different for the nn, pp, or np cases. 

Here, we will first describe our DBHF calculation of asymmetric matter (Section II). 
In Section III we will present and discuss results for the EOS at various levels of asym- 
metry, which confirm the findings from Ref. [18]. We then compare high-density pressure 
calculations to most recent experimental observations. In the context of non-relativistic 
calculations, we also consider a popular potential based on chiral perturbation theory. In 
Section IV, we explore the sensitivity of the G-matrix at various positive energies to increas- 
ing values of the neutron excess parameter. Conclusions and plans for future applications 
are presented in Section V. 



II. DESCRIPTION OF THE CALCULATION 

Asymmetric nuclear matter can be characterized by the neutron density, p n , and the 
proton density, p p . It is also convenient to define the total density p = p n + p p and the 
asymmetry (or neutron excess) parameter a = Pn ~ pp , Clearly, a=0 corresponds to symmetric 
matter, and a=l to neutron matter. 

In terms of a and the average Fermi momentum, related to the total density in the usual 
way, 

the neutron and proton Fermi momenta can be expressed as 
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kp = k F (l + a) 1/s 



(2) 



and 

k p F = k F (l - af 3 , (3) 

respectively. 

We use the Thompson relativistic three-dimensional reduction of the Bethe-Salpeter 
equation. The Thompson equation is applied to nuclear matter in strict analogy to free- 
space scattering and reads, in the nuclear matter rest frame, 

9ij(q, Q, P, (4)o) = v ij(<f, 0) 

f d 3 K -*m*m* Q tJ (K,P) 

+ J (a^V* ' K) ^f (4)o - 4( P, k) 9 ^ * P ' ( ^ )o) 

(4) 

where ij—nn, pp, or np, and the asterix signifies that medium effects are applied to those 
quantities. In Eq. (4), q, (f, and K are the initial, final, and intermediate relative momenta, 



and E* = y (m*) 2 + K 2 . The momenta of the two interacting particles in the nuclear matter 
rest frame have been expressed in terms of their relative momentum and the center-of-mass 
momentum, P, through 

pJ-±±h (5) 

K = ^ (6) 

The energy of the two-particle system is 

e* J (P,K) = e*(P,K) + e*(P,K) (7) 

and (e*j)o is the starting energy. The single-particle energy e* includes kinetic energy and 
potential energy, see below. The Pauli operator, Q, prevents scattering to occupied states. 
To eliminate the angular dependence from the kernel of Eq.(4), it is customary to replace 
the exact Pauli operator with its angle-average. Detailed expressions for the Pauli operator 
in the case of two different Fermi momenta are given in Appendix A. It is also customary to 
introduce an average center-of-mass momentum [5]. Definitions and details are also provided 
in Appendix A. 



With the definitions 



and 



_ m* m 3 



V* = m t v*. ^ (9) 
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one can rewrite Eq.(4) as 



(10) 

which is formally identical to the non-relativistic equation. We make use of the definitions 
Eqs.(8-9) throughout this paper, which is why all formulas involving the opertor G appear 
identical to their non-relativistic equivalent. 

The goal is to determine self-consistently the nuclear matter single-particle potential 
which, in our case, will be different for neutrons and protons. To facilitate the description of 
the numerical procedure, we will use a schematic notation for the neutron/proton potential 
(while the corresponding detailed expressions are reported in Appendix B). We write, for 
neutrons, 

U n = U np + U nn (11) 

and for protons 

Up = Up n + Upp (12) 

where each of the four pieces on the right-hand-side of Eqs. (11-12) depends on the appro- 
priate G-matrix (nn, pp, or np) from Eq.(10). Clearly, the two equations above are coupled 
through the np component and so they must be solved simultaneously. Furthermore, the 
G-matrix equation and Eqs. (11- 12) are coupled through the single-particle energy (which 
includes the single-particle potential). So we have three coupled equations to be solved self- 
consistently. As done in the symmetric case [9], we parametrize the single particle potential 
for protons and neutrons (Eqs. (11-12)) in terms of two constants, Us,% and Uy,i, (the scalar 
and vector potential) through 

nrn* 

U i (k i ) = -^U s , i + U v , i . (13) 

For the purpose of facilitating the connection to the usual non-relativistic framework [5], it 
is customary to redefine Us,i and Uy,i in terms of two other constants defined as 

m* = rrii + U s>i (14) 

and 

U ,i = U s , + U Vtl . (15) 

The subscript "i" signifies that these parameters are different for protons and neutrons. 
Starting from some initial values of m* and Uo^, the G-matrix equation is solved and a first 
approximation for Ui(ki) is then obtained. This solution is again parametrized in terms of 
a new set of constants, and the procedure is repeated until convergence is reached. The 
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effective masses for neutrons and protons, m*, obtained through this procedure are shown 
in Figs. (1) as a function of the Fermi momentum at various levels of asymmetry. 
Finally, the energy per neutron or proton in nuclear matter is calculated from 

Ei =< Ti > + <Ui> . (16) 

(See Appendix C for details.) The EOS, or energy per nucleon as a function of density, is 
then written as 

_/ \ Priori + PpCp /irr\ 

e{pn,pp) = (17) 

or 

-d \ (l + a)e„ + (1 -a)e p 

e(k F ,a) = . (18) 



III. PROPERTIES OF THE EQUATION OF STATE 

The NN potential used in this work is the relativistic OBEP from Ref. [20] which uses 
the Thompson equation and the pseudo-vector coupling for the n and r] mesons. 

The EOS as obtained from our DBHF calculation is displayed in Fig. 2 (upper panel), as 
a function of kp and for values of a between and 1. The symmetric matter EOS saturates 
at k F ps 1.4/m -1 with a value of 16.7 MeV, in good agreement with the empirical values. 
For the compression modulus of saturated symmetric matter, defined as 

K = ^W^-r (19) 

we find a value of 233 MeV. This is in excellent agreement with the recent empirical deter- 
mination of 225±15 MeV [21]. 

As the neutron density increases (the total density remaining constant), the EOS becomes 
increasingly repulsive and the minimum shifts towards lower densities. As the system moves 
towards neutron matter, the "energy well" gets more and more shallow, until, for a larger 
than about 0.8, the system is no longer bound. 

We also show for comparison the EOS based on the BHF calculation, see Fig. 2, lower 
panel. Also in this case we have performed a self-consistent calculation based on the same 
realistic force. However, no medium modifications are included in the potential to account 
for the proper Dirac structure of the nucleons in nuclear matter. As a consequence of that, 
saturation of symmetric matter is obtained at a much higher density, a well-known problem 
with the conventional approach. 

In Fig. 3, we take a different look at the EOS. There, we plot the quantity e(k F , a) — 
e(kp, 0) versus a 2 . Clearly, the behaviour is linear, that is 

e(kp , a) — e(k F , 0) = e s a 2 , (20) 

or parabolic versus a. This linear behaviour, shared with the non-relativistic predictions, 
see lower panel of Fig. 3, is reminiscent of the asymmetry term in the familiar semi-empirical 
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mass formula. By reading the slope of each line, we can then predict the nuclear symmetry- 
energy, defined as 

_ l d 2 e(k F ,a) 

e s - - d 2 a l"=o- ( 21 ) 

This is shown in Fig. 4, where the solid curve is the prediction from the DBHF model 
and the dashed corresponds to the BHF calculation. The DBHF prediction at saturation 
density is about 30 MeV. The behaviour of the symmetry energy is most controversial at 
high densities, where even the trend depends strongly on the type of interaction used [22]. 
Energetic reactions induced by heavy neutron-rich nuclei have been proposed as a mean to 
obtain crucial information on the high-density behaviour of the nuclear symmetry energy 
and thus the EOS of dense neutron-rich matter [22]. Analyses of collision dynamics have 
been done at the National Superconducting Cyclotron Laboratory to extract EOS-sensitive 
observables such as the elliptic flow [19]. In this way, empirical constraint has been obtained 
for the pressure of symmetric and neutron matter. 

The pressure of neutron-rich matter is defined in terms of the energy/particle as 

P M= ^J^l (22) 

(We notice that Eq. (20) facilitates the calculation of the pressure for any a.) 

In Fig. 5, we show the pressure obtained from the DBHF-based EOS for symmetric 
matter and neutron matter. The empirical constraints have been provided by the authors 
of Ref. [19]. The area within the dashed lines indicates EOS consistent with the elliptic 
flow measurements reported in that work. For the case of neutron matter, two different 
parametrizations for the symmetry energy are assumed in the analysis. The two pressure 
contours in the bottom panel of Fig. 5 correspond to the weakest (lower contour) and the 
strongest (higher contour) density dependence for e s (p) proposed by Prakash et al. [23]. 

It is fair to say that DBHF predictions produce a reasonable amount of repulsion. In 
contrast, relativistic mean field theories (RMFT) tend to generate too much pressure [19], 
as can be expected from the higher compression modulus usually predicted by RMFT. 
While very repulsive EOS such as obtained from RMFT as well as weakly repulsive EOS 
(compression modulus of 167 MeV) are ruled out by the analysis, EOS which become softer 
at densities higher than three times saturation density (possibly due to a phase transition 
to quark matter) are not excluded. 



Non-relativistic calculations with a high-precision chiral potential 

The predictions shown above are typical of meson-exchange potentials with moderate 
strength in the tensor force (as indicated by a deuteron D-state probability of about 5%). 
In-medium predictions, however, can be very sensitive to the nature of the chosen two- 
body force, and can be so even if the interactions under consideration predict very similar 
free-space observables. 

Recently, a new generation of NN potentials based on chiral perturbation theory (%PT) 
has begun to emerge [24,25]. These interactions may become increasingly popular in the 
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future as a way to describe low-energy processes in terms of effective degrees of freedom 
while respecting the symmetries of QCD. Thus we felt it would be interesting to consider 
one of these potentials as well in calculations of symmetric and asymmetric matter. In 
this subsection, we show basic EOS properties for the chiral potential of Ref. [24], which is 
quantitative in its description of NN data up to 300 MeV. 

Chiral potentials are non-relativistic in nature, and thus unsuitable for a DBHF calcula- 
tion. For comparison, we include predictions from a modern high-precision meson-exchange 
potential (the CD-Bonn from Ref. [26]). (The symmetry energy for various high-precision 
potentials has been calculated in Ref. [27].) CD-Bonn, which uses pseudoscalar coupling 
for the pion and has been developed within the framework of the non-relativistic Lippman- 
Schwinger equation, will also be used, as it is appropriate, in a conventional Brueckner 
calculation. 

When evaluationg these predictions, one must keep in mind that consideration of many- 
body forces is a crucial issue in non-relativistic calculations (whereas some three-body forces 
are effectively included in the DBHF framework [28]). The chiral potential, in particular, 
requires the inclusion of three-body forces for consistency within its order. However, we are 
interested in a baseline comparison between the two interactions as it will reveal differences, 
if any, originating at the two-body level. In fact, the free-space predictions of the two 
potentials under our consideration are of very similar quality. However, the momentum 
structure of chiral potentials can be quite different than the one of conventional meson- 
exchange potentials, with the chiral interactions being "softer" due to cutoffs applied to 
eliminate high-momentum components, as required by the nature of the chiral expansion. 
Such differences are likely to show up in G-matrix calculations. 

EOS properties from the chiral potential in comparison with CD-Bonn are displayed in 
Figs. 6-8. Generally, the behaviour is similar to the one of a low tensor force potential, 
such as Bonn A, when used in a non-relativistic context (the chiral potential predicts a 
deuteron D-state probability of 4.93%). Although similar at moderate densities, the two 
sets of predictions start differing considerably for higher values of the Fermi momentum. 
The best overview is provided by the symmetry energy, Fig. 8, where we see the two curves 
being very close to each other up to about saturation density, after which the chiral curve 
seems to get a bit steeper, but than crosses over to become softer at the higher densities. 

Due to the cutoff applied to the chiral potential which suppresses momentum components 
higher than 460 Mev/c [24], it's legitimate to ask up to which value of the Fermi momentum 
its predictions can be deemed reliable. As a guideline for this estimate, we will take the 
case of symmetric matter, which is obtained from the formulas in the appendices setting the 
proton and neutron Fermi momentum equal to each other. Evaluation of the energy/particle 
requires integration over the relative momentum (denoted as Kq in Appendices B-C) up to 
&f, which is therefore the maximum value of Kq. With the cutoff of 460 Mev/c being 
equivalent to 2.3 fm -1 , the chiral potential predictions are realistic for the values of kp 
shown in Figs. 6-8. 

IV. IN-MEDIUM INTERACTIONS 

The G-matrix calculated from Eq. (10) has the dependence (for a given energy and 
momentum) , 
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Gf = G^\k F ,a), (23) 

where f3 is an appropriate set of quantum numbers, and ij refers to the type of nucleons. 
We now wish to solve Eq. (10) for scattering states. Envisioning a physical situation where 
protons are the incident particles, we then have pp and pn possibilities, which are calculated 
separately with the appropriate effective masses and Pauli blocking operator. 

As is well known since the birth of nuclear physics, the pp and pn interactions in the 
same state would be identical except for charge dependence. However, in the asymmetric 
medium, these interactions are physically different due to the different Fermi momenta of 
protons and neutrons. In Fig. 9, the neutron and proton Fermi momenta from Eqs. (2-3) 
are shown versus the asymmetry parameter a for fixed (average) Fermi momentum. 

To which extent the asymmetry in the proton/neutron ratio in a nucleus will affect 
proton-nucleus scattering is equivalent to the question: how sensitive is the interaction to 
the asymmetry degree of freedom? We can certainly get insight into this issue from the 
(infinite nuclear matter) G-matrix, which, upon proper mapping into coordinate space and 
application of the local density approximation, will yield the effective interaction for finite 
nuclei [29,30]. When "folded" over the nucleus, the interaction yields the optical potential 
used to describe elastic scattering, or, in the case of inelastic scattering, it becomes the 
transition potential connecting the initial state to the target excited states. 

At present, scattering calculations on halo nuclei are being done with phenomenological 
interactions where the isospin asymmetry is accounted for through different parametrizations 
for pp and pn [31]. The Melbourne group [32] has produced a very comprehensive set of (non- 
relativistic) proton-nucleus scattering calculations, including some exotic nuclei. However, 
no attempt is made to account for the asymmetry in the interaction. Clearly, additional 
calculations would be helpful. 

For this sensitivity study, we have selected one S-wave and one P-wave with T=l, so 
that pp and pn cases are both existing and can be compared. The 1 S'o wave will give us 
insight into the central force, while the 3 Pi wave contains strong sensitivity to the spin-orbit 
component. (The central and spin-orbit pieces of the effective interaction play the most 
important role in proton-nucleus elastic scattering and inelastic scattering to natural-parity 
transitions.) 

In Figs. 10-11, we show the real and imaginary parts of the (on-shell) 1 5'o matrix element 
at three different energies. In each case, the quantity is plotted versus the (average) Fermi 
momentum, for three values of a. The same legend applies to Fig. 12-13 for the 3 Pi matrix 
element. 

In general, the (smaller) imaginary part is much more sensitive to the degree of asym- 
metry. The presence of an imaginary part is determined by the occurrence of a singularity 
in the scattering equation, Eq. (10). That occurrence will depend on the kinematics and 
the value of the Fermi momenta. We see for instance that at 50 MeV the imaginary part 
vanishes at higher densities for larger values of a. This is particularly striking in the pp 
case, due to the drastic lowering of the proton Fermi momentum as seen from Fig. 9. 

It is interesting to notice how the a-dependence shows up differently in the pp/pn cases. 
The dependence on a comes in through the p/n effective masses (see Fig. 1) and of course 
Pauli blocking, which acts differently on pp or pn states (see Appendix A for details). Pauli 
blocking in pp scattering is controlled by just the proton Fermi momentum, whereas pn 
scattering is controlled by both the neutron and proton Fermi momenta, which vary as seen 
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from Fig. 9. As a result, the scattering probability in the pn states exhibits overall a weaker 
dependence on a. 

We would expect a similar trend for in-medium pp or pn cross sections, a very useful 
concept for transport model calculations [33], where both the mean field and nucleon-nucleon 
collisions play an important role. For instance, we would expect in-medium pp cross sections 
to be more sensitive than pn cross sections to increasing asymmetry. By the same token, 
nn cross sections would be strongly Pauli-blocked in neutron-rich matter. We will look into 
this aspect in a future work. 

Finally, as the energy increases, the details of Pauli blocking become much less important, 
and so does the a-dependence, although the overall density dependence remains strong. (At 
higher energies, relativistic density-dependent effects, which are included, become the most 
important.) Experiments at low energy would be best to reveal sensitivity of the interaction 
to the asymmetry between neutron and proton distributions. 

V. SUMMARY AND CONCLUSIONS 

We have presented a microscopic calculation of the equation of state of nuclear matter 
when protons and neutrons have different Fermi momenta. The calculation is self-consistent 
and parameter-free, in the sense that no parameter of the NN force is adjusted in the 
medium. 

As expected, the single-particle energy moves up to less attractive values to merge with 
the neutron matter equation of state when the proton density approaches zero. The depen- 
dence of the EOS on the neutron excess parameter is clearly linear as a function of a 2 . We 
make predictions for the nuclear symmetry energy and observe a large discrepancy between 
the relativistic and the non-relativistic predictions at high density. Our fundings are in 
agreement with those from Ref. [18]. 

A relatively simple way to relate the microscopic EOS directly to structural properties of 
finite nuclei is the use of a mass formula [14], where the "volume" term is directly related to 
the EOS. In that context, it has been pointed out [14] that, unlike proton densities, neutron 
densities are very sensitive to the EOS model. At the same time, precise data on neutron 
radii and neutron skins are not readily available to discriminate clearly among models. Even 
for a nucleus such as 208 Pb, for which a fairly large database exists, determinations of the 
neutron skins differ considerably from model to model [34]. We are presently in the process 
of using our EOS in calculations of neutron radii and neutron skins. 

The relative simplicity of a homogeneous infinite system makes nuclear matter calcula- 
tions a convenient starting point for the determination of an effective interaction suitable for 
finite nuclei. Together with the local density approximation, this approach has been used 
extensively and with success for proton-nucleus scattering. With that in mind, we have 
taken a closer look at the scattering matrix in asymmetric matter. 

Ultimately, the goal is to study nuclei with high levels of asymmetry, about which very 
little is known. Coherent effort from both the experimental and the theoretical side is 
necessary in order to combine reliable models of the density dependence of the effective 
nuclear force with reliable structure information. Hopefully, the RIA facility will be available 
in the near future to answer open questions and help us understand the physics of the weakly 
bound systems that are expected to exist at the limits of the nuclear chart. 
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APPENDIX A: ELIMINATION OF THE ANGULAR DEPENDENCE 
1. Angle-averaged Pauli operator 

The definition of the Pauli operator is, for the case of two different Fermi momenta 



Q= { 





kn < kp \ 


if < 


or 


< 


k u < kp J 




k/j, > kp I 


1 if < 


and / 




k u > kp J 



(fi^u = 1,2) 



(Al) 



In the angle-average approximation, one replaces the exact Q operator with its average 
over all angles for fixed P and K. That is, one defines 



Q=<Q>- 



J Q(K, P, 9)dQ _ 1 r 
7^ ~ 2 hi 



r Q(P, K, 6)d9 

J6i 



(A2) 



The variables P, K, k^ are as defined previously, 9 is the angle between P and K, and 
Q ^ for 9 1 < 9 < 2 . 

Several cases can occur, depending on the values of K , P , k F and k F . Those are: 



Q = 



ro k 2 < \{{k F f + (k F ) 2 ) - p 2 \ 

1 (P-K) 2 >(k F ) 2 

^{{K + P) 2 -{kp) 2 ) (kp) 2 < (P - K) 2 < (k F ) 2 > 
_i_ (K* + p2 _ i^ k ny + (a .p )2 ^ otherwise 



(A3) 



2. Averaged center- of- mass momentum 

To simplify the integration leading to the single-particle spectrum and the energy per 
particle (see Appendices B and C), we introduce the average center-of-mass momentum. 
This is defined as the root-mean-square value of P for two particles with the constraint 
that their relative momentum is K (case a), and, in addition, that one of the particles has 
momentum k^ (case b). 
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a. K fixed 

Definition of average cm. momentum (K fixed): 

jf F dh jf F dk 2 P 2 5(K h\) 



p z 



J^dh J kF dk 2 5(K - \ \k x -k 2 \ 



To simplify the final expressions, we introduce the following notation 
x = k n F + K y = k p F -K s = k F -K t = kp + K 
The final expression is then 



p2 
r AV 



l{k F f + K 2 K<k p F 
§K(s 5 + y 5 ) + ±(ty + sx) 3 + §(s 6 + y 6 ) - (ty 5 + sx 5 ) 



= < _ _ 

\K(s 3 + y 3 ) + \{ty + sx) 2 + (s 4 + y A ) - 2{ty 3 + 

.0 



otherwise 



{k F + k p F ) < 2K 



b. K and fixed 
Definition of average cm. momentum (K and k^ fixed): 

So F dk 2v 5{K — | \kft~ k v \) 



kfj, < kp 



Piv = < 



l(3(k,) 2 + (kp) 2 -4Kk,) 

n 



ku > kp 



2^K <C kp k^ 
kp k^ < IK < kp -\- k^ 
otherwise 



p2 _ 
r AV — 



\(3(k,) 2 + (k F ) 2 -AKk,) 



kji kp < IK < kp -\- k^ 
otherwise 



where ji = 1, 2 and i=n or p. 
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APPENDIX B: NUCLEAR MATTER SINGLE-PARTICLE POTENTIAL 



The single-particle potential in nuclear matter, Ui (i=n or p), is defined in the usual way 



Ui(k^) = < I'J'j G lJ t',Vj - Vj^i >= 

j<A 



< ViVj 1 \G in \n i v h - v h iii > + < iiiV h \G %v \[LiV h - > 

ji<N j 2 <Z 

The explicit evaluation of this expression in partial-wave decomposition yields 



(Bl) 



8A 



M L,a 



rku+Ko 

K dK / PdPGr L n (P;K ,K ) 



k"+k 



+ K ° dK ° -, v , r -,. 



2 "(J 



PdPGr L n (P;K ,K ) 



/ n < — > p 



ftp 



kp 



(a stands for the quantum numbers J, S and T). 

Using the averaged- momentum approximation (see Appendix A), 

G a L i j (P; K , K ) « G a L i j (P A V, K , K ), 
allows to simplify the integral over P, with the final result being 

8A 



[ ' 2 K dK G a L i n (PAv; K , K ) [ k " +K ° PdP 

JO ' J\k„-K \ 



(B2) 



(B3) 



fc2 _i_(un \2 



I* K dK Gr L n (pAv,K ,K ) r * F ~ K °PdP 

l\k^-K \ 



+ 



n < — > p 
k F < — ► k F 



16A 



7T 



L,a 



K 2 dK Gr L n (pAV,K ,K ) 



2k„, J V 



(k n F ) 2 ~ kj 



K (K -k,) Gr L n (pAV,K ,K ) 



+ 



n < — > p 
— kl 



(B4) 



The coefficients contain the isospin dependence and are equal to 



/t-T=0,1 _ -t-T=0,1 _ 1 . q~0 _ n-0 _ n . 

2 np 2 pn 2 ' nn pp ' 

r 1 = r 1 = i 

nn pp 
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APPENDIX C: ENERGY PER NUCLEON IN NUCLEAR MATTER 



The energy per neutron/proton is: 



2iV j<Ai<N 



(CI) 



1 



22 j<Ai<Z 

which, in view of Eq.(Bl), can be written as 



(C2) 



e n =< T n > +— Y, ^(M 

ZiV i<N 



(C3) 



e p =< T p >+—J2 Ui{k„) 



i<Z 



(C4) 



By applying again the partial-wave decomposition and the average momentum approxi- 
mation 



G£?'(P; K , K ) « G a L i j (PAV, K , K ) 
the energy per neutron/proton is finally 



(C5) 



Ei =<Ti> + 



6A 



f / 2 K" dK Gt P \P AV ;K ,K 



fc F + fc F 



{kpY + (£^) 3 ((&j£) 2 — (&f) 2 )_ ^ K| _ "^ ((/ c ^,) 2 + (A;^) 2 ) 



3 2 



3 jo 



2-3-^+ 3 



(C6) 



where, again, i = n or p. 
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FIG. 1. Upper panel: neutron effective mass as a function of the Fermi momentum and for 
increasing values of a between and 1. The upper curve corresponds to a = 1. Lower panel: 
proton effective mass. The upper curve corresponds to a = 0. 
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FIG. 2. Upper panel: energy per nucleon as a function of the Fermi momentum at differente 
values of the asymmetry parameter (in steps of 0.1) from symmetric matter (lowest curve) to 
neutron matter (highest curve). The predictions are obtained from DBHF calculations. Lower 
panel: corresponding predictions from BHF calculations. 
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FIG. 3. Upper panel: the left-hand side of Eq. (20) versus a 2 (DBHF model) for increasing values 
of the Fermi momentum from 0.9 fm^ 1 (lowest curve) to 1.7 fm^ 1 . Lower panel: corresponding 
predictions from BHF calculations. 
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FIG. 4. Nuclear symmetry energy as a function of the Fermi momentum. The solid line is 
the prediction from DBHF calculations, while the dashed line is obtained with the conventional 
Brueckner approach. 
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FIG. 5. Upper panel: pressure (in MeV//m 3 ) in symmetric matter as a function of density (in 
units of saturation density). The area within the dashed lines indicates pressure values consistent 
with the data. Lower panel: pressure in neutron matter. The two contours correspond to two 
parametrizations of the asymmetry term used in the analysis, as explained in the text. The 
predictions are from DBHF calculations. 
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FIG. 6. Upper panel: energy per nucleon as a function of the Fermi momentum at differente 
values of the asymmetry parameter (in steps of 0.2) from symmetric matter (lowest curve) to 
neutron matter (highest curve). The predictions are obtained from BHF calculations with the 
chiral potential. Lower panel: corresponding predictions with CD-Bonn. 
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FIG. 7. Upper panel: the left-hand side of Eq. (20) versus a 2 (obtained with the chiral potential) 
for increasing values of the Fermi momentum from 0.9 fm~ l (lowest curve) to 2.2 /m _1 . Lower 
panel: corresponding predictions with CD-Bonn. 
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FIG. 8. Nuclear symmetry energy as a function of the Fermi momentum. The solid line is the 
prediction obtained with the chiral potential, while the dashed line is obtained with CD-Bonn. 
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FIG. 9. The proton (lower line) and neutron (upper line) Fermi momenta as a function of the 
asymmetry parameter for fixed value of the (average) Fermi momentum (solid line) . 
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FIG. 10. Real part of the np and pp Sq matrix element as a function of the Fermi momentum 
at different energies and for different levels of asymmetry. 
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FIG. 11. Imaginary part of the np and pp So matrix element as a function of the Fermi mo- 
mentum at different energies and for different levels of asymmetry. 
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FIG. 12. As in Fig. 10, for the 3 Pi matrix element. 
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FIG. 13. As in Fig. 11, for the 3 Pi matrix element. 
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